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IMPULSIVE  LOADING  OF  A  SIMPLY  SUPPORTED  CIRCULAR  PLATE 


by 

Norman  Jones* 


Abstract 

It  is  clear  from  a  survey  of  literature  on  the  dynamic  deformation  of 
rigid-plastic  plates  that  most  work  has  been  focused  on  plates  in  which 
either  membrane  forces  or  bending  moments  alone  are  considered  important , 

t 

while  the  combined  effect  of  membrane  forces  and  bending  moments  on  the  be¬ 
havior  of  plates  under  static  loads  is  fairly  well  established.  This  arti¬ 
cle  ,  therefore ,  is  concerned  with  the  behavior  of  circular  plates  loaded  dy¬ 
namically  and  with  deflections  in  the  range  where  both  bending  moments  and 
membrane  forces  are  important. 

A  general  theoretical  procedure  is  developed  from  the  equations  for 
large  deflections  of  plates  and  a  simplified  yield  condition  due  to  Hodge. 

It  can  be  shown  that  these  general  equations  may  be  reduced  to  give  the  pre¬ 
dictions  of  Onat  and  Haythornthwaite  for  static  loading  and  to  those  of  Wang 
for  the  dynamic  case  with  bending  moments  only  if  the  rigid,  perfectly  plas¬ 
tic  material  from  which  the  plate  is  made  yields  according  to  Tresca  and  any 
strain  rate  effects  are  disregarded. 

The  results  obtained  when  solving  the  governing  equations  for  the  par¬ 
ticular  case  of  a  simply  supported  circular  plate  loaded  with  a  uniform  im¬ 
pulsive  velocity  are  found  to  compare  favorably  with  the  corresponding  ex¬ 
perimental  values  recorded  by  Florence.  Although  this  study  is  of  general 
interest,  it  is  thought  in  particular  that  it  should  assist  in  the  interpre¬ 
tation  of  the  dynamic  biaxial  stress-strain  characteristics  of  materials  re¬ 
corded  on  diaphragms  fitted  in  impact  tubes. 
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defined  by  equations  (52),  (53)  and  (55) 

arbitrary  constants  of  integration 

defined  by  equations  (45)  and  (46) 

plate  thickness 

impulse  per  unit  area  of  plate 

Bessel  function  of  order  n 

defined  by  equation  (54) 

a  H2/4 
o 

radial  and  circumferential  bending  moments  per 
unit  length 

a  H 
o 

radial  and  circumferential  membrane  forces  per 
unit  length 

transverse  shear  force  per  unit  length  of  plate 
outside  radius  of  plate 
principal  radii  of  curvature 

initial  velocity  of  plate 

Neumann's  Bessel  function  of  the  second  kind  of 
zero  order 

uniform  distributed  pressure  per  unit  area  of 
undeformed  plate 

dimensionless  bending  moments  M  /M  ,  M„/M 

&  r  o’  0  o 

dimensionless  membrane  forces  N  /N  ,  N  /N 

r  o’  0  o 


P 


-k  sin<f> 


Notation  (continued) 


q  -k  cos<}> 

r  radial  coordinate  of  plate 

t  time 
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yV  R2 
o 

12M 

o 

duration  of  first  stage  of  deformation 

duration  of  entire  deformation  corresponding 
to  the  greatest  value  from  equation  (48) 

displacement  in  direction  r  of  undeformed 
plate 

transverse  deflection  perpendicular  to  unde¬ 
formed  plate 

final  deformation 

r/R 

a  r/R 
s 

1  +  E 

r 

s  roots  of  equation  J  (a  )  =  0 

o  s 

u  +  r 

yV  2R/M 
o  o 

W“>1/2 

H/R 

radial  and  circumferential  strains 

circumferential  coordinate  lying  in  plate 

a  /R 
s 

mass  per  unit  area  of  plate 


Notation  (continued) 


P 

a 

o 

$ 


radius  of  hinge 

yield  stress  in  simple  tension 


slope  of  the  mid-plane  of  a  plate  measured  in  a 
plane  which  passes  through  r  =  0  and  is  per¬ 
pendicular  to  the  plate  surface 


Introduction 


During  the  second  world  war  Taylor  (1),  Hudson  (2),  Richardson  and  Kirk¬ 
wood  (3)  and  others  (4)  were  the  first  to  conduct  experimental  and  theoretical 
studies  into  the  influence  of  dynamic  loads  on  the  behavior  of  thin  disks  or 
circular  plates.  Further  work  on  membranes  has  been  reported  by  Frederick  (5), 
who  proposed  a  mechanism  of  behavior  somewhat  similar  to  that  of  Hudson  (2), 
and  Griffith  and  Vanzant  (6),  who  recorded  dynamic  load-carrying  capacities 
significantly  greater  than  the  corresponding  static  values.  It  is  also  evident 
from  their  results  that,  at  high  loading  rates,  elements  of  such  plates  tend  to 
move  in  a  transverse  sense  which  ensues  in  smaller  circumferential  strains  than 
those  expected  from  static  tests. 

Hopkins  and  Prager  (7)  contributed  some  general  information  on  the  dynamic 
behavior  of  plates  and  studied  the  particular  problem  of  a  simply  supported 
circular  plate  subjected  to  a  uniformly  distributed  load  which  is  brought  on 
suddenly,  and,  after  being  kept  constant  during  a  certain  interval  of  time,  is 
removed  suddenly.  The  plate  is  made  of  a  rigid,  perfectly  plastic  material, 
which  is  assumed  to  obey  the  Tresca  yield  condition  and  associated  flow  rule. 
Any  membrane  forces  which  may  arise  during  deformation  are  disregarded.  Making 
the  same  assumptions  and  using  a  similar  but  more  complicated  mechanism  of  be¬ 
havior,  Florence  (8)  solved  the  problem  of  a  circular  plate  clamped  around  its 
outer  edge  and  loaded  with  a  central  rectangular  pulse.  Wang  and  Hopkins  (9) 
studied  the  behavior  of  a  circular  plate  with  a  transverse  velocity  imparted 
instantaneously  to  the  entire  plate  except  at  the  built-in  outer  edge,  where 
the  velocity  is  zero.  Wang  (10)  found  that  this  analysis  simplifies  consider¬ 
ably  when  the  plate  is  simply  supported  around  its  outer  edge.  Florence  (11) 
subjected  some  simply  supported  circular  plates  to  uniformly  distributed 


0 
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impulses  and  found,  particularly  for  large  impulses,  that  the  theoretical 
analysis  in  (10)  overestimated  considerably  the  final  deformed  shapes.  This 
discrepancy,  however,  arises  mainly  because  the  nature  of  the  assumptions  in¬ 
volved  in  the  development  of  the  analysis  in  (10)  limits  application  of  the 
results  to  plates  with  small  final  deformations. 

Perzyna  (12)  examined  the  influence  of  pulses  of  arbitrary  shape  by  de¬ 
veloping  further  the  theory  of  Hopkins  and  Prager  (7)  to  show  that,  for  a 
given  impulse,  the  character  of  the  pressure-time  function  has  little  influ¬ 
ence  on  the  final  shape  of  the  plate.  In  passing,  however,  it  is  worth  not¬ 
ing  that  Sankaranarayanan  (13)  showed  that  the  pulse  shape  influenced  con¬ 
siderably  the  final  deformation  of  plastic  spherical  caps  subjected  to  impact 
pressure  loads,  while  Symonds  (14)  found  somewhat  less  sensitivity  to  the 
pulse  shape  in  beams. 

Hopkins  (15)  developed  a  more  general  theory  for  plates  loaded  trans¬ 
versely  with  non -symmetrical  loads,  but  disregarded  any  membrane  forces  and 
solved  no  particular  problems. 

Shapiro  (16),  who  was  the  first  to  examine  the  dynamic  behavior  of  an¬ 
nular  plates,  studied  the  problem  of  a  circular  plate  supported  rigidly 
around  an  inner  radius  and  loaded  with  a  circular  ring  of  impulsive  velocity 
around  its  outer  edge.  Recently,  Florence  (17)  reviewed  this  problem  in 
order  to  assess  the  relative  contributions  of  membrane  forces  and  bending 
moments  to  the  formation  of  the  final  deformed  shape.  It  was  found  that  a 
solution  which  considered  interaction  between  the  circumferential  membrane 
force  and  circumferential  bending  moment  was  closer  to  the  experimental 
values  than  the  solutions  for  bending  moment  only  and  membrane  force  only, 
both  of  which  overestimated  considerably  the  final  deformations.  Witmer, 
Balmer,  Leech  and  Pian  (18)  developed  a  numerical  method  and  a  computer 
program,  the  predictions  of  which  compare  favorably  with  experimental 
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values  recorded  for  the  large  dynamic  deformations  of  beams,  rings,  plates 
and  shells. 

Boyd  (19)  reconsidered  the  problem  of  dynamic  deformation  of  a  circular 
membrane  and  solved  the  governing  equations  numerically  for  a  general  form  of 
symmetrical  pressure  loading  which  permits  variations  across  the  plate  and 
with  time.  Although  any  contributions  arising  from  bending  moments  were  neg¬ 
lected,  this  method  predicted  results  similar  to  the  corresponding  values 
computed  in  Ref.  (18).  Boyd  (19)  and  Frederick  (5)  investigated  the  deforma¬ 
tions  of  membranes  made  from  a  strain-hardening  material  and  discovered  that 
a  simplified  rigid,  perfectly  plastic  analysis  provides  a  remarkably  accurate 
approximation  to  the  true  behavior.  Munday  and  Newitt  (20)  examined  carefully 
the  dynamic  behavior  of  clamped  copper  membranes,  and  presented  some  photographs 
taken  from  a  high-speed  cine  film  recorded  during  the  passage  of  a  hinge  which 
formed  around  the  outer  edge  of  a  plate  at  the  first  instant  of  impact  and 
travelled  inwards  towards  the  center,  where  it  remained  until  the  plate  came  to 
rest.  For  an  impulsive  loading  which  imparts  to  the  plate  a  velocity  increas¬ 
ing  from  zero  at  the  first  instant  of  impact  to  a  maximum  value  at  the  center 
some  time  later,  Munday  and  Newitt  observed  that  the  final  deformed  shapes  are 
almost  identical  for  all  plates  irrespective  of  thickness,  diameter  and  load. 
Johnson,  Poynton,  Singh  and  Travis  (21)  performed  some  underwater  explosive 
stretch-forming  experiments  on  clamped  circular  blanks,  and  measured  the  thick¬ 
ness  strains  across  the  plates  and  final  deformed  shapes  for  various  charge 
strengths . 

It  is  of  passing  interest  to  note  that  the  method  of  Martin  and  Symonds 
(22)  can  be  used  to  predict  the  maximum  deflections  and  time  bounds  for  rigid- 
plastic  plates  loaded  impulsively.  Although  the  time  bounds  compare  rather 
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well  with  those  predicted  in  (9,  10),  the  deflections  are  overestimated  by 
one  third  for  the  simply  supported  case  and  somewhat  less  for  the  plate  with 
clamped  edges.  It  is  worth  emphasizing,  however,  that  these  results  are  ob¬ 
tained  from  a  few  lines  of  arithmetic  whereas  the  solution  in  (9)  in  particu¬ 
lar  is  time  consuming. 

Cooper  and  Shifrin  (23)  and  Haythorthwaite  and  Onat  (24)  measured  the 
static  load -carrying  capacity  of  initially  flat  circular  plates, and  observed 
that  the  bending  only  solution  of  Hopkins  and  Prager  (25)  underestimates  con¬ 
siderably  the  load  which  could  be  supported  if  deflections  of  the  order  of  the 
plate  thickness  or  larger  are  permitted.  In  order  to  explain  the  strengthen¬ 
ing  effect  under  static  loads,  Onat  and  Haythomthwaite  (26)  considered  both 
bending  moments  and  membrane  forces  in  an  analysis  which  reduces  to  the  bend¬ 
ing  only  theory  of  Ref.  (25)  at  very  small  deflections  and  to  the  behavior  as 
a  membrane  at  large  deflections.  Symonds  and  Mentel  (27)  have  shown  that  a 
simple  analysis  of  plastic  beams  constrained  axially  and  loaded  with  a  trans¬ 
verse  pressure  impulse  becomes  unrealistic  if  deflections  of  the  order  of  the 
beam  thickness  are  permitted  and  the  influence  of  axial  forces  is  disregarded. 
The  authors  describe  the  transition  of  clamped  and  simply  supported  beams  from 
an  initial  simple  behavior  to  a  final  stage  in  which  deformations  are  governed 
primarily  by  catenary  effects. 

It  is  clear  from  the  foregoing  survey  of  literature  that  most  effort  has 
been  concentrated  on  the  dynamic  deformation  of  plates  in  which  either  mem¬ 
brane  forces  (1-5,  19)  or  bending  moments  (7-10,  12,  15,  16)  alone  are  be¬ 
lieved  to  be  important.  Moreover,  with  the  exception  of  the  computing  work 
in  (18)  and  the  analysis  by  Florence  (17)  for  an  annular  plate,  no  investiga¬ 
tions  have  been  conducted  into  the  interaction  effects  between  membrane  forces 
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and  bending  moments,  while  for  static  loading  the  behavior  during  all  stages 
is  fairly  well  understood  (26). 

The  object  of  this  present  article,  therefore,  is  an  attempt  to  link  the 
two  distinct  stages  of  plastic  strain  and  describe  the  behavior  of  plates  dy¬ 
namically  loaded  with  deflections  in  the  range  where  both  bending  moments  and 
membrane  forces  are  important.  Subsequently,  it  is  found  that  the  theoretical 
analysis  presented  herein  predicts  with  reasonable  accuracy  the  final  deforma¬ 
tions  recorded  by  Florence  (11)  on  a  simply  supported  circular  plate  subjected 
to  a  uniform  impulse. 


Equilibrium  Equations 

Introducing  the  inertia  terms  into  the  equilibrium  equations  developed 
by  Reissner  (28)  for  the  finite  deflection  of  a  circular  plate  loaded  stati¬ 
cally  gives 

(aeV'  “  VNe  "  Ve  Q/Rr  +  Vep  "  pae™'  “  wVe'fi  =  0  (1) 

(aeQ)'  +  VecVRr  +  W  +  arV  +  vVe'”  '  y(Vw'  =  0  (2) 

SV’  '  VMe  “  Ve  Q  =  0  (3) 

provided  the  rotary  inertia  effect  is  disregarded,  and 

a  =  1  +  e  , 
r  r 

aQ  =  r  +  u  =  r(l  +  e0)  , 

1/Rr  =  ♦'/(!  +  er)  , 
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1/Rg  =  sin<j>/r  , 

( * )  =  |r-  (  )  ,  t  =  time 


(  )’  = 


y 


The  positive  directions  of  the  various  quantities  are  indicated  in  Fig.  1. 

Since  we  are  interested  in  the  interaction  between  membrane  forces  and 
bending  moments,  then  we  may  follow  Reissner  (28)  and  limit  equations  (1)  - 
(3)  to  plates  having  small  strains  and  deflections  which  are  not  too  large. 
Thus,  using 


a.  =  r,  a  =1,  1/R  =  1/R.  =  sin<|>/r,  and  a  '  =  cos<J> 

o  r  r  o  d 

equations  (1)  -  (3)  become 


(rN  )'  -  cos(J>N.  -  rffi'Q  +  rp  -  yrww'  -  yr  cos<|>u  =  0  (4) 

r  0 

(rQ)’  +  rN  <f>'  +  sinifrN  +  rq  +  yr  cos<J>w  -  yruw'  =0  (5) 

r  o 

(rM  )'  -  cos<f>Ma  -  rQ  =  0  (6) 

r  o 

If  it  is  further  assumed  that  cos<|>  =  1  and  sin<|>  =  -w'j,  and  recognized 
that  r<f>'Q  and  rNr«{i'  are  small  in  comparison  with  the  remaining  terms,  then 
equations  (4)  -  (6)  can  be  rewritten  in  the  form 


m  '  +  n  -  n„  =  -rkw’/R  +  yrww'/N  +  yrii/N  (7) 

r  r  0  o  o  o 

rm  "  +  2m  ’  -  mQ'  -  4n0w'/H  =  rk/M  -  yrw/M  +  yriiw'/M  (8) 

r  r  0  0  o  o  o 


Ve  =  Ve/No  . 


m  =  M  ,„/M  , 

r  0  r’0  o  ’ 


where 
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and  the  second  order  terms  <j>'w’  are  disregarded.  In  view  of  the  observa¬ 
tions  of  Griffith  and  Vanzant  (6)  that  the  material  of  a  circular  plate  tends 
to  move  in  a  transverse  sense  at  high  rates  of  dynamic  loading,  then  we  as¬ 
sume  hereafter  that 

u  =  0  (9) 


Yield  Condition 

It  is  necessary  to  employ  a  four-dimensional  yield  surface  between  the 
generalized  stresses  n^,  n^,  m^  and  m^  ,  in  order  to  solve  the  equilibrium 
equations  (7)  and  (8)  for  a  plate  made  from  a  rigid-plastic  material.  Onat 
and  Prager  (29)  have  derived  equations  for  the  load -carrying  capacity  of 
shells  of  revolution  made  from  a  rigid,  perfectly  plastic  material  which 
obeys  Treses' s  yield  condition  and  the  associated  flow  rule.  However,  these 
relations  are  very  difficult  to  handle  for  other  than  simple  problems.  Hodge 
(30),  therefore,  proposed  a  two-moment  limited  interaction  surface  which 
maintains  all  interaction  between  force  and  force  and  moment  and  moment ,  but 
neglects  all  interactions  between  force  and  moment .  If  it  is  assumed  that 
yielding  is  controlled  by  the  Tresca  criteria,  the  result  is  a  linear  surface 
in  four-dimensional  space;  and  the  equations  of  the  twelve  planes  are  indica¬ 
ted  in  Fig.  2. 

The  yield  condition  proposed  by  Hodge  (30)  was  used  to  solve  the  equilib¬ 
rium  equations  (7)  and  (8)  for  a  simply  supported  plate  loaded  with  a  static 
pressure  distributed  uniformly  over  a  central  circular  area.  This  method  pre¬ 
dicted  a  variation  of  load  carrying  capacity  very  similar  to  that  forecasted 
by  Onat  and  Haythorthwaite  (26)  for  deflections  greater  than  H/2. 
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Mechanism  of  Deformation 


The  notion  of  a  travelling  hinge  has  been  used  successfully  for  estab¬ 
lishing  the  dynamic  behavior  of  beams  and  cantilevers  (27,  31,  etc.). 
Furthermore,  Munday  and  Newitt  (20)  observed  that  at  the  first  instant  of 
impact  a  hinge  develops  at  the  supports  of  a  circular  plate  and  travels  in¬ 
ward  towards  the  center, where  it  remains  until  the  plates  come- to  rest.  It 
seems  reasonable,  therefore,  to  employ  herein  the  similar  mechanism  of  be¬ 
havior  proposed  by  Wang  (10)  for  a  simply  supported  circular  plate  loaded 
impulsively. 

Thus  at  t  $  0  ,  the  plate  travels  with  an  initial  uniform  velocity  v  , 

o 

and  a  hinge  is  assumed  to  form  at  the  support  radius  r  =  R  at  the  instant 
t  =  0  .  Some  time  later  (t),  the  hinge  will  reach  a  radius  p(t)  ,  where 
0  <  p(t)  <  R  ,  and  divide  the  plate  into  two  distinct  zones  as  indicated  in 
Fig .  3  where 


w  =  V  ,  for  0  $  r  *  p(t) 


(10) 


and 


w  = 


V  (R  -  r)/(R 
o 


p) 


for  p(t)  $  r  S  R 


(11) 


It  is  assumed  that  the  kinetic  energy  of  the  plate  for  the  moment  when 
the  travelling  hinge  reaches  radius  p  =  0  is  dissipated  during  a  second  stage 
of  deformation  in  which  the  hinge  remains  stationary  at  r  =  0  . 


Expression  for  Slope  (w1) 

If  p(t)  $  r  $  R  ,  then  using  (10)  and  (11)  the  transverse  deflection  of 


the  plate  at  time  t  is 


(12) 


w  =  V  t(r)  + 
o 


3w 

f  dt 


t(r) 


where  is  the  initial  uniform  velocity  of  the  plate  at  t  =  0  ;  t(r)  is 

the  time  taken  for  a  travelling  hinge  to  reach  radius  r  ,  provided  p(t)  $ 
r  S  R  ,  while  t  is  the  time  at  which  the  hinge  reaches  radius  p  . 

Letting 

—  -  R  T 
St  "  R  T 

where  R  is  a  function  of  r  alone ,  and  T  is  a  function  of  t  alone  per¬ 
mits  (12)  to  be  rewritten 


w  =  V  t(r)  +  R 
o 


T  dt 


t(r) 


whence , 


!?■».{?<«»»♦!?(*> 


T  dt  +  R  |- 
Sr 


t(r) 


f 


•t  ^ 

T  dt 

t(r) 


(13) 


Analysis  of  a  Simply  Supported  Circular  Plate  Loaded  Impulsively 
(i)  First  Stage  of  Deformation 


It  can  be  shown  that  the  yield  conditions , 


m  =  mA  =  -1  ,  n  =  n.  =  1  ,  for  0  $  r  $  p 
r  o  r  o 


(14) 


and 


mQ  =  -1  ,  0  s  m  $  -1  ,  n  =1,  0  <  nQ  $  1  ,  for  p  $  r  $  R  (15) 

o  r  r  o 
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indicated  in  Fig.  2,  are  consistent  with  equations  (10)  and  (11).  Equations 
(10)  and  (14)  satisfy  exactly  equations  (7),  (8),  (9)  and  (13),  while  (7), 
(8),  (9)  and  (15)  give 


provided 


4r  3w 
H  3r 


32w  kr2 

.  .  2  M 
3t^  o 


p(t)  $  r  $  R  ,  and  (w' )2  <<  1  . 


(16) 


If  3w/3r  =  0  and  k  =  0  ,  then  for  the  first  stage  of  deformation  it  can 
be  shown  that  equations  (11)  and  (16)  yield  the  following  time  function. 


t  =  tlci  -  £  -  (|)2  ♦  (|)3] 


(17) 


which  is  identical  to  that  obtained  by  Wang  (10)  who  neglected  the  effect  of 
membrane  forces. 

Since  one  would  expect  the  influence  of  the  bending  moments  to  dominate 
over  the  action  of  the  membrane  forces  during  the  first  stage,  then  the  form 
of  the  time  function  for  3w/3r  +  0  will  be  taken  as 


t  .  tci  -  |  -  (£)2  +  (|)3]  (18) 

where  t  is  a  constant  to  be  determined  later.  Differentiating  equation  (18) 
with  respect  to  time  gives 

1/(R  -  p)  =  -pt( 3p  +  R)/R3  (19) 

which  allows  (13)  to  be  rewritten  in  the  form 

~  =  -2-  (^§—  +  Rp  -  - - Rr)  ,  for  p(t)  S  r  s  R  .  (20) 

3r  ^3  2  z 
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Substituting  equations  (11)  and  (20)  into  (16)  yields 


4V  t 

—  (r2  m  ' )  =  — —  (^-  p2r  +  Rpr  -  r3  -  Rr2) 
ar  r  r3h  2  2 


y^V  r2(R  -  r) 
Mq(R  -  p)2 


(21) 


which  when  integrated  with  the  condition  that  m  =  -1  at  r  =  p  gives 


4  V  t 


m 


o  /  3 


R3H 


/3  j  .  R  r3  Rr2  ,  3p4  .  Rp3  ,  Rp 

(it  P  r  +  o  pr  -  Q - —  +  8F  +  6F  “  p  "  ? 


ypv, 


£  rB.  r,2  II  +  £p_  PI  i 

2  v6  "  12  3r  "  4r  ~  2  3' 


MQ(R-p) 


(22) 


Since  =  0  at  r  =  R  ,  equation  (22)  yields 


7V  t  yV  R2 
o  o 

6H~  “  Iff 
o 


+  1  =  0,  at  p  =  0 


(23) 


the  solution  of  which  can  be  written  as 

2  1/2 


V  t  =  R 
o 


(.  M  +  I 

1  7R  2 


2yHV 


18HM 


(- 


7M 


-)  (1  + 


7yV  2R2 
o 


1/2 

'  i 


(24) 


or. 


3  oF  1/2  q  p 

4  *-!'♦<&> 


(25) 


where , 


e  =  H/R 


B  =  yV  2R/M 
M  o  o 


(26) 

(27) 


and 


A  =  Vt/R 


(28) 
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V  t  is  the  deflection  of  the  center  of  the  plate  at  the  conclusion  of 
the  first  stage  (p  =  0). 

It  can  be  shown  that  the  deflected  form  of  the  plate  at  t  =  t  is 

”  ■  V[1  '  y  <s>’  -  y  ^ 


(ii)  Second  Stage  of  Deformation 


Now, one  would  expect  the  membrane  forces  to  dominate  over  the  bending 
moments  during  the  second  stage  of  deformation  throughout  which  the  hinge 
remains  stationary  at  r  =  0  .  Thus#  assuming  =  Q  =  0  and  consider¬ 

ing  the  case  when  k  =  0  ,  equations  (1),  (2)  and  (9)  become 


(r  n  )  -  cos<|>nfl  =  yrww'/N 

r  o  o 


r^'n  +  sin<Jm„  =  -vircos<f>w/N 
r  o  o 


Utilizing  the  same  yield  condition  as  used  previously  for  p  $  r  $  R 
(15)  permits  equations  (30)  and  (31)  to  be  rewritten  as  one  equation. 


w"  +  w'/r  =  Pw/N 


since , 


(rn  ) '  =  m  '  +  cos<bn 
r  r  r 

2 

cos$  =  1  ,  sin$  =  -w'  ,  and  (w')  <<  1 


The  general  solution  of  (32)  may  be  written  in  the  following  form  after 


Bowman  ( 32 ) 
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w  =  {A  J  (X  r)  +  B  Y  (X  r)}  {C  cos( 
os  os 


n  In" 

—  X  t)  +  D  sin  (  —  X  t)} 

y  s  4  y  s 


(33) 


for  t  {  t  j  ts 


or 


pr  jw~ 

w  =  {E  cos  (  —  X  t)  +  F  sin  (  —  X  t 

s  y  s  s  y  s 


)}J  (X  r) 
o  s 


(34) 


in  order  to  maintain  finite  deflections  at  r  =  0  . 

The  deflection  w  must  be  zero  at  the  simply  supported  edge  for  all 
values  of  t  .  This  restriction  on  (34)  yields  the  relation, 


Xs  =  R  ’  8  *  1»  2  »  ’ 


(35) 


where  a  are  the  roots  of  the  equation  J  (a  )  =  0  .  Thus, 
s  ^  o  s  ’ 


w  =  l  {E  cos(y  t )  +  F  sin(y  t)}  J  (X  r) 

IS  s  s  s  OS 


(36) 


where 


.  Po 
Ys  J  y 


(37) 


The  deflection  w  and  velocity  w  given  by  (36)  and  its  derivative 
must  match  the  values  acquired  by  the  plate  at  the  end  of  the  first  stage 
given  by  (29)  and  (11)  with  p  =  0  and  t  =  t  . 

Hence , 

2  3  00 

V  t[l  -  ^ - =  l  {E  cos(y  t )  +  F  sin(y  t)}  J  (a  x)  (38) 

O  4  4  -  s  s  s  s  OS 

S  =  1 


and 
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V  (1  -  x)  =  l  {-Y  E  sin(Y.t)  +  Y  F  cos(y  t)]  J  (a  x)  (39) 

O  S  S  S  SS  SJOS 

S=1 

where , 

x  =  r/R  .  (40) 


Multiplying  both  sides  of  (38)  and  (39)  by  x  J  (a  x)  dx  and  inte- 

o  s 

grating  over  the  limits  x  =  1  and  x  =  0  ,  it  can  be  shown  with  the  aid 
of  Lommel’s  integrals  (Ref.  32),  that 

“1  2  3 

x  V  t(l  -  - - ^-)  J  (o  x)dx  (41) 

O  2  2  o  S 

o 

•1 

x  V  (1  -  x)  J  (a  x)dx  (42) 

o  os 

o 

where 


H  = 
s 


J,Z(a  ) 


G  = 
s 


J,2(a  ) 


G  =  E  cos(y  t)  +  F  sin(y  t) 
s  s  s  s  s 


(43) 


and 


H  =  -YE  sin(Y_t)  +  y^F  cos(Yot ) 
s  s  s  s  s  s  s 


(44) 


If  y  =  a  x  ,  equations  (41)  and  (42)  become 
s 


V  o 

G  =  - 2 -  {i3  j  (a  )  -  — 

s  -2 \  1  s  „  2 


J,  a_ 

1  s  s 


s 


JQ(y)dy  ) 


(45) 


and 


2V 


H  = 
s 


Jl2(as)as3  J 


JQ(y)dy 


(46) 
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where  values  of  the  Bessel  functions  are  tabulated  in  references  (33)  and 
(34).  Schmidt  (35)  tabulates  the  integral 

f«s 

JQ(y)dy  for  10  $  y  $  50  . 

■  o 

The  plate  will  come  to  rest  when  w  =  0  .  Hence , 


l  l- 


s=l 


Y  E 
s  s 


sin(y  t  )  +  v  F  cos(y  t  )}  J  (X  r)  =  0 
ss  ss  ss  OS 


tan(Ysts)  =  Fs/Es  . 

Thus  the  final  deformed  shape  w^.  at  t  =  tg  is 

"f  ■  i,  k  ♦  ri  jo<v> 

s=l 

or 

*f*  l  fl  ♦  (Hs/rs)2  J0(Xsr) 
s=l 

which  may  be  expressed  in  the  form 


»f  *  R  I  AS  i(4Bs>2  ♦  K  Jo<Xsr> 
S=1 

where 

As  =  1/{J12(0S)  o|}  , 

Bs  =  13  J^a.)  -  9  Ks/«!  , 


(47) 


(48) 


(49) 


(50) 


(51) 


(52) 


(53) 
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K 


s 


s  JQ(y)dy 
'o 


(54) 


and 


C 

s 


K  /a 

S  5 


(55) 


Conclusions 

A  general  theoretical  procedure  which  retains  both  membrane  forces  and 
bending  moments,  but  neglects  strain  rate  effects,  is  developed  herein  in 
order  to  describe  the  behavior  of  rigid,  perfectly  plastic  circular  plates 
loaded  with  axisymmetrical  dynamic  loads.  It  can  be  shown  that  this  analy¬ 
sis  may  be  reduced  to  the  predictions  of  Onat  and  Haythorthwaite  (26)  for 
static  loading  of  simply  supported  plates  and  to  those  of  Wang  (10)  for  the 
dynamic  case  with  bending  moments  only. 

It  is  evident  from  Fig.  4  that  the  maximum  deflection  at  r  =  0  predic¬ 
ted  by  the  first  five  terms  of  equation  (51)  compares  favorably  over  a  wide 
range  of  uniformly  distributed  impulses  with  the  experimental  values  ob¬ 
served  by  Florence  (11)  on  simply  supported  circular  aluminium  plates  having 
R/H  =  16  .  Theoretical  predictions  of  the  maximum  plate  deflections  are 
presented  in  dimensionless  form  in  Fig.  5  for  various  plate  parameters  and 
magnitudes  of  impulse . 

It  is  thought  that  the  method  developed  could  be  used  to  describe  the 
behavior  of  plates  having  other  support  conditions  and  various  dynamic  loads. 
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FIG.  2  YIELD  CONDITION  AFTER  HODGE  (30) 
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oo  EXPERIMENTAL  RESULTS  FOR  AL606I-T6 
R/H  =16  (AFTER  FLORENCE  (II)) 

- BENDING  ONLY  THEORY  (WANG  (10)) 

- MEMBRANE  -  BENDING  THEORY 

FIG.  4  COMPARISON  OF  THEORY  WITH 
EXPERIMENTAL  RESULTS  FOR  AN 
IMPULSIVELY  LOADED  SIMPLY 
SUPPORTED  PLATE 
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FIG. 5  DEFLECTION  -  IMPULSE  RELATION  FOR  SIMPLY 
SUPPORTED  CIRCULAR  PLATES 


